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ABSTRACT 

Models  of  the  very  early  Universe,  including  inflationary  models,  are  argued  to  produce  varying 
universe  domains  with  different  values  of  fundamental  constants  and  cosmic  parameters.  Using 
the  cosmological  hydrodynamical  simulation  code  from  the  eagle  collaboration,  we  investigate 
the  effect  of  the  cosmological  constant  on  the  formation  of  galaxies  and  stars.  We  simulate 
universes  with  values  of  the  cosmological  constant  ranging  from  A  =  0  to  Aq  x  300,  where 
Ao  is  the  value  of  the  cosmological  constant  in  our  Universe.  Because  the  global  star  formation 
rate  in  our  Universe  peaks  at  f  =  3.5  Gyr,  before  the  onset  of  accelerating  expansion,  increases 
in  A  of  even  an  order  of  magnitude  have  only  a  small  effect  on  the  star  formation  history 
and  efficiency  of  the  universe.  We  use  our  simulations  to  predict  the  observed  value  of  the 
cosmological  constant,  given  a  measure  of  the  multiverse.  Whether  the  cosmological  constant 
is  successfully  predicted  depends  crucially  on  the  measure.  The  impact  of  the  cosmological 
constant  on  the  formation  of  structure  in  the  universe  does  not  seem  to  be  a  sharp  enough 
function  of  A  to  explain  its  observed  value  alone. 

Key  words:  cosmological  parameters  -  dark  energy  -  inflation. 


1  INTRODUCTION 

Cosmological  inflation,  it  has  been  argued,  naturally  predicts  a  vast 
ensemble  of  varying  universe  domains,^  each  with  different  cosmic 
conditions  and  even  different  fundamental  constants  (see  the  review 
of  Linde  2017).  A  typical  mechanism  for  generating  these  universes 
is  as  follows  (Guth  2007).  The  inflation  fleld  undergoes  quantum 
fluctuations,  and  so  we  might  expect  some  parts  of  the  universe 
to  still  be  inflating  while  other  parts  have  entered  a  post-reheating 
‘big  bang’  phase.  The  universe  as  a  whole  consists  of  post  big  bang 
universes  fllled  with  ordinary  matter  and  radiation,  surrounded  by 
an  ever-inflating  background. 

In  evaluating  such  models,  predicting  what  we  would  expect 
to  observe  is  necessarily  tied  to  where  observers  are  formed  in  the 
multiverse.  In  this  instance,  anthropic  reasoning  is  inevitable  (Carter 
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^  For  simplicity,  we  call  such  regions  ‘universes’. 


1974;  Carr  &  Rees  1979;  Davies  1983;  Barrow  &Tipler  1986).  With 
different  cosmic  and  fundamental  constants  in  different  parts  of  the 
multiverse,  the  values  we  expect  to  observe  are  unavoidably  tied  to 
their  ability  to  support  the  complexity  required  by  life. 

These  multiverse  models  could  successfully  explain  the  flne  tun¬ 
ing  of  the  universe  for  life:  small  changes  in  their  values  can  suppress 
or  erase  the  complexity  upon  which  physical  life  as  we  know  it,  or 
can  imagine  it,  depends.  The  scientiflc  literature  on  the  flne  tuning 
of  the  universe  for  life  has  been  reviewed  in  Hogan  (2000),  Barnes 
(2012),  Schellekens  (2013),  MeiBner  (2015),  and  Lewis  &  Barnes 
(2016).  For  example,  as  pointed  out  by  Davies  &  Unwin  (1981), 
Sakharov  (1984),  Linde  (1984),  Banks  (1985),  Linde  (1987),  Wein¬ 
berg  (1987),  and  Weinberg  (1989),  only  a  small  subset  of  values 
of  the  cosmological  constant  (A)  permit  structure  to  form  in  the 
universe  at  all.  Universes  in  which  the  cosmological  constant  is 
large  and  positive  will  expand  so  rapidly  that  gravitational  struc¬ 
tures,  such  as  galaxies,  are  unable  to  form.  Large  negative  values 
will  cause  space  to  recollapse  rapidly,  also  preventing  the  formation 
of  galaxies. 

If  inflation  creates  a  huge  number  of  variegated  universe  domains, 
then  a  structure-permitting  value  of  the  cosmological  constant  will 
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probably  turn  up  somewhere.  Any  observers  will  see  a  universe 
with  at  least  some  structure.  In  this  way,  the  seemingly  improbable 
suitability  of  our  universe  for  life  is  rendered  more  probable. 

As  Weinberg  (1987)  noted,  we  can  test  a  particular  multiverse 
model  via  its  prediction  of  the  distribution  of  universe  properties. 
Observers  will  inhabit  universes  drawn  in  a  highly  biased  way  from 
the  population  of  universes,  but  we  can  calculate  the  typical  prop¬ 
erties  of  a  universe  that  contains  observers.  In  this  way,  we  can 
calculate  the  likelihood  of  our  observations,  and  so  compare  mul¬ 
tiverse  models.  For  example,  a  model  in  which  99  per  cent  of  ob¬ 
servers  measure  a  value  of  the  cosmological  constant  as  large  as  our 
value  should  (other  things  being  equal)  be  preferred  over  a  model 
in  which  only  I  per  cent  of  observers  make  such  a  measurement. 
Whether  these  consistency  tests  can  give  absolute  (rather  than  just 
relative)  support  to  the  idea  of  a  multiverse  is  the  subject  of  some 
debate  (Ellis  &  Silk  2014;  Barnes  2017). 

To  test  the  relative  merits  of  multiverse  models  in  this  way,  we 
need  to  know  how  life,  or  at  least  the  cosmic  structures  that  are  the 
likely  preconditions  for  life,  depends  on  the  fundamental  constants 
of  nature  and  cosmic  parameters.  In  the  case  of  the  cosmological 
constant,  the  large-scale  structure  of  the  universe  is  most  directly 
affected.  Galaxies  are  the  sites  of  star  formation,  and  stars  provide 
both  a  steady  source  of  energy  and  the  heavier  elements  from  which 
planets  and  life  forms  are  made. 

Within  an  anthropic  approach,  we  can  also  shed  light  on  the 
coincidence  problem:  we  live  at  a  time  in  the  Universe  when  the 
energy  density  of  the  cosmological  constant  and  the  energy  den¬ 
sity  of  matter  are  within  a  factor  of  2  of  each  other  (Lineweaver 
&  Egan  2007).  The  coincidence  problem  has  motivated  a  search 
to  alternative  modification  to  gravity  that  might  explain  the  value 
of  the  cosmological  constant  more  naturally.  Although,  alternative 
models,  such  as  quintessence  can  explain  why  the  relative  densi¬ 
ties  of  matter  and  cosmological  constant  densities  track  each  other, 
fine  tuning  of  the  model  parameters  is  still  required  to  explain  their 
observed  similarity  (Zlatev  &  Steinhardt  1999;  Zlatev,  Wang  & 
Steinhardt  1999;  Dodelson,  Kaplinghat  &  Stewart  2000;  Chimento 
et  al.  2003). 

Investigations  of  the  effect  of  the  cosmological  constant  on  galaxy 
formation  have  thus  far  relied  on  analytic  models  of  increasing 
levels  of  sophistication.  Efstathiou  (1995)  located  galaxies  at  the 
peaks  of  the  smoothed  density  field  of  the  universe,  and  found  that 
-  assuming  that  the  cosmological  constant  is  positive,  observers 
should  expect  to  see  ~  0.67-0.9.  Peacock  (2007)  extended  this 
approach  to  negative  values  of  the  cosmological  constant,  finding  a 
significant  probability  that  <  0  is  observed.  These  approaches 
have  been  extended  by  Garriga  &  Vilenkin  (2000),  Garriga,  Livio 
&  Vilenkin  (2000),  Tegmark  et  al.  (2006),  Bousso  &  Leichenauer 
(2009),  Bousso  &  Leichenauer  (2010),  Piran  et  al.  (2016),  Sudoh 
et  al.  (2017),  and  Adams  et  al.  (2017). 

The  modem  approach  to  galaxy  formation  uses  supercomputer 
simulations  that  incorporate  the  effects  of  gravity,  gas  pressure,  gas 
cooling,  star  formation,  black  hole  (BH)  formation,  and  various 
kinds  of  feedback  from  stars  and  BH  accretion.  It  has  been  long 
known  that  feedback  is  very  important  to  explaining  the  star  for¬ 
mation  history  of  our  Universe;  models  without  feedback  are  too 
effective  at  forming  stars,  compared  to  observations  (White  &  Rees 
1978;  Dekel  &  Silk  1986;  White  &  Erenk  1991;  Somerville  &  Dave 
2015).  One  of  the  key  ingredients  that  has  allowed  this  progress  is 
the  inclusion  of  realistic  models  for  the  impact  of  feedback  from  the 
growth  of  BHs.  All  successful  models  now  demonstrate  the  need  for 
active  galactic  nuclei  (AGNs)  as  an  additional  source  of  feedback 
that  suppresses  the  formation  of  stars  in  high-mass  haloes  (Benson 


et  al.  2003;  Bower  et  al.  2006;  Croton  et  al.  2006).  Although  this 
idea  was  initially  developed  using  semi-analytic  models,  this  has 
now  been  confirmed  in  a  wide  range  of  numerical  simulations  (e.g. 
Dubois  et  al.  2016;  Bower  et  al.  2017;  Pillepich  et  al.  2018). 

Here,  we  will  use  the  eagle  project’s  galaxy  formation  code  to 
calculate  the  effect  of  the  cosmological  constant  on  the  formation 
of  structure  in  different  post-inflation  universes.  Each  of  our  models 
will  be  practically  indistinguishable  at  early  times,  including  nucle¬ 
osynthesis  and  the  epoch  of  recombination.  Their  histories  diverge 
at  later  times  due  to  the  onset  of  cosmological  constant-powered 
accelerating  expansion.  In  Section  2,  we  describe  the  eagle  galaxy 
formation  code  and  the  suite  of  simulations  that  we  have  run.  In  Sec¬ 
tion  3,  we  describe  the  effect  of  changing  the  cosmological  constant 
on  the  global  accretion  and  star-forming  properties  of  the  universe. 
Section  4  looks  at  the  effect  on  an  individual  galaxy,  and  its  relation 
to  its  environment.  In  Section  5,  we  use  our  simulations  to  derive 
prediction  from  models  of  the  multiverse. 

2  GALAXY  FORMATION  SIMULATION  CODE 

The  Virgo  Consortium’s  eagle  project  {Evolution  and  Assembly 
of  GaLaxies  and  their  Environment)  is  a  suite  of  hydrodynamical 
simulations  that  follow  the  formation  of  galaxies  and  supermassive 
BHs  in  cosmologically  representative  volumes  of  a  standard  A 
cold  dark  matter  (ACDM)  universe.  The  details  of  the  code,  and 
particularly  the  sub-grid  models,  are  described  in  Schaye  et  al. 
(2015),  and  are  based  on  the  models  developed  for  OWLS  (Schaye 
et  al.  2010),  and  used  also  in  GIMIC  (Crain  et  al.  2009)  and  cosmo- 
OWLS  (Le  Brun  et  al.  2014).  The  simulations  code  models  the 
effect  of  radiative  cooling  for  11  elements,  star  formation,  stellar 
mass-loss,  energy  feedback  from  star  formation,  gas  accretion  on 
to  and  mergers  of  supermassive  BHs,  and  AGN  feedback. 

The  initial  conditions  for  the  eagle  simulations  were  set  up  us¬ 
ing  a  transfer  function  generated  using  cams  (Lewis,  Challinor  & 
Lasenby  2000)  and  a  power-law  primordial  power  spectrum  with 
index  n^  =  0.9611.  Particles  were  arranged  in  a  glass-like  initial 
configuration  were  displaced  according  to  second-order  Lagrangian 
perturbation  theory  (Jenkins  2010). 

BHs  are  seeded  in  all  dark  matter  haloes  with  masses  greater 
than  10^^  h~^  Mq  =  1.48  x  10^^  Mq.  The  halo  finding  algorithm 
is  described  in  Schaye  et  al.  (2015);  in  short,  the  code  regularly  runs 
the  friends-of-friends  (EoE)  finder  (Davis  et  al.  1985)  with  link¬ 
ing  length  0.2  on  the  dark  matter  distribution.  When  analysing  the 
simulations  in  following  sections,  we  are  interested  in  membership 
with  any  halo,  rather  than  distinguishing  substructures,  so  we  use 
the  EOE  algorithm  to  identify  haloes. 

2.1  Cosmological  parameters  and  scale  factor 

We  need  to  choose  the  cosmological  parameters  for  our  simula¬ 
tion.  The  problem  with  the  standard  set  of  cosmological  parameters 
(^^m,  ^a,  h)  is  that  they  are  all  time  dependent.  In  the  model 
universes  that  we  will  consider,  there  is  no  unique  ‘today’  at  which 
we  can  compare  sets  of  parameters.  We  follow  Tegmark  et  al.  (2006) 
by  defining  cosmological  parameters  that  are  constant  in  time.  We 
use  only  one  time-dependent  parameter,  which  is  cosmic  time  t.  The 
constant  parameters  are  listed  in  Table  1.  Note  that  the  cosmological 
constant  (A)  and  its  associated  energy  density  are  related  linearly, 
A  =  %iiG  . 

How  do  we  solve  the  Eriedmann  equations,  given  the  dimen¬ 
sionless  cosmological  parameters  in  Table  1,  so  that  we  can  derive 
the  usual  cosmological  parameters  for  the  simulation?  We  have  the 
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Table  1.  Free  parameters  in  the  FLRW  model,  defined  so  that  they  are  constant  in  time, 
at  least  since  very  early  times.  The  measured  value  derives  from  the  Planck  Collaboration 
XVI  (2014)  cosmological  parameters,  as  used  by  the  eagle  project:  a%,ns,Y)  = 

(0.307,  0.693,  0.048  25,  0.6777,  0.8288,  0.9611,  0.248). 


Parameter 

Measured  value 

Pa 

Cosmological  constant  energy  (mass)  density 

5.98  X  10“^^  kg  m“^ 

Baryon  mass  per  photon  po/ny 

1.01  X  10“^^  kg  m“^ 

Cold  dark  matter  mass  per  photon  Pc/ny 

5.43  X  10“^^  kg  m“^ 

K 

Dimensionless  spatial  curvature  (in  Planck  units)  k/a^TQ 

\k\  <  10-60  ^0 

1  1 

freedom  to  choose  ‘today’,  i.e.  we  can  rescale  a{t)  to  make  aito)  =  1 
for  any  time  ^o-  A  useful  way  to  proceed  initially  is  to  define  to  be 
the  time  at  which  the  energy  densities  of  the  cosmological  constant 
and  matter  are  equal.  Then,  we  calculate  the  matter  densities, 

Pm,o  =  Pk  (1) 


^ 


^b,o  -  ^Pm,o  ; 

Sm 


Pc,o  —  -^Pm,o- 

Sm 


Then,  we  calculate  the  photon  number  density  at  to,  and  from  it  the 
CMB  temperature  (Tq)  and  the  radiation  (photons  and  neutrinos) 
energy  density. 


—  Pb,o/?b  —  Pc,o/?c  —  Pm,o/fm 


^K,0  — 


2^(3)  rksTo 
71^  V  he 


Pr,o  =  gi^^BTo 


ksTo 

he 


where 


^  =  2  +  2 


7 

8 


3 


(4) 

(5) 

(6) 


We  can  then  solve  the  Friedmann  equation. 


if  =  - 


1  da 


SnG 


a  dt 


(Pm  ~\~  Pr  ~\~  P A  ~\~  Pk) 


(7) 


SttG 


(Pm,o^  +  Pr,o<^  Pa-\-  PkfiCi  ), 


(8) 


3  f  ksTo 

where  Pk  =  —k -  - 

^  SttG  V  ^ 

We  can  calculate  the  critical  density,  Pent, 0  =  3 //q /SttG  =  yOm,o  + 

Py,o  +  Pa  +  Pk,o  and  then  the  usual  cosmological  parameters 
=  Pm,o/Pcrit,  0.  and  similarly  for  12c,  ^a,  and  12k.  With 

these  parameters,  FLRW  codes  can  solve  the  Friedmann  equations.^ 
Having  solved  the  Friedmann  equations  for  a(t),  we  can  rescale 
to  change  the  time  of  ‘today’  to  be  any  other  time  (tg):  a^ewit)  = 
a(t)/a(tQ),  and  recalculate  the  various  density  parameters  appropri¬ 
ately.  We  will  describe  our  choices  for  the  normalization  of  a(t)  in 
Sections  2.3  and  2.4. 


2.2  Initial  conditions  and  sub-grid  physics 

We  use  the  same  initial  conditions  for  each  simulation.  For  the  range 
of  cosmological  constants  we  consider  here,  there  has  been  minimal 
effect  on  the  evolution  of  the  universe  at  the  start  of  the  simulation. 
Specifically,  we  use  the  same  initial  conditions  for  the  SPH  particles 
in  physieal  eoordinates:  in  the  eagle  code,  like  its  gadget  ancestor, 
we  need  to  convert  code  quantities  into  physical  quantities  taking 
into  account  the  initial  scale  factor  (a^)  and  the  Hubble  parameter 
(h)  of  the  original  simulation:  distance  (<iphys  =  ciih~^dcode)^  velocity 

(b’phys  —  and  mass  (wip^ys  —  h  ^code)- 

We  must  also  be  careful  regarding  parameters  in  our  sub-grid 
physics  recipes.  The  sub-grid  physics  of  the  eagle  code  has  been 
checked,  and  the  necessary  parameters  rescaled  as  necessary  to  keep 
the  same  physieal  values.  We  also  discovered  a  few  cases  in  which 
it  was  assumed  that  pa  ^  0,  which  needed  to  be  remedied  for  the 
test  runs  below. 

Note  the  assumptions  that  we  are  making  when  we  change  the 
cosmological  constant,  but  keep  the  physical  parameters  of  the  sub¬ 
grid  model  unchanged.  This  is  potentially  worrisome,  given  that 
these  parameters  are  often  inferred,  not  from  first  principles,  but 
by  calibrating  against  observations  of  galaxy  populations  in  our 
Universe.  Our  assumptions  are  twofold.  First,  we  assume  that  the 
sub-grid  model  is  sufficiently  sophisticated  that  it  captures  the  rel¬ 
evant  physics.  For  example,  we  assume  that  star  formation  in  any 
cosmology  occurs  when  the  local  density  is  sufficiently  high.  It  is 
appropriate  to  apply  such  a  model  to  other  universes.  Secondly, 
we  assume  that  the  parameters  inferred  from  observations  are  the 
same  as  would  be  inferred  from  a  first-principles  calculation;  they 
do  not  depend  on  the  cosmological  constant  for  such  small-scale 
processes.  For  example,  the  loeal  matter  density  above  which  star 
formation  occurs  should  only  depend  on  conditions  within  10-100 
pc  scale  molecular  clouds,  far  below  cosmological  scales.  We  can 
plausibly  use  the  same  threshold  for  different  cosmologies. 

Using  the  same  sub-grid  parameters  would  create  a  problem  only 
if  our  overall  cosmology  is  wrong,  for  it  could  be  the  case  that  we 
have  inferred  the  wrong  value  of  some  sub-grid  parameter  to  par¬ 
tially  compensate  for  an  incorrect  expansion  history  of  the  universe. 
In  this  case,  of  course,  the  entire  eagle  simulation  suite  would  need 
to  be  redone,  as  would  almost  every  other  cosmological  simulation. 
We  will  leave  that  worry  for  another  day. 


^  There  are  two  potential  complications.  If  we  consider  a  universe  with  no 
cosmological  constant  (pa  =  0),  then  the  choice  of  ‘initiaT  matter  density  is 
effectively  arbitrary.  Secondly,  if  the  universe  recollapses,  then  it  may  never 
reach  the  time  at  which  pm,o  =  Pa-  The  most  general  way  to  find  some 
matter  density  at  which  we  can  apply  the  technique  above  is  to  write 
the  Friedmann  equation  in  terms  of  the  CMB  temperature  Tfi.  We  can  then 
solve  for  the  CMB  temperature  at  turnaround  H{t)  =  0,  and  from  this  calcu¬ 
late  the  minimum  matter  density  of  the  universe. 


2.3  Testing  our  modifications 

The  freedom  to  choose  ‘today’  to  in  our  simulation  gives  us  a  way 
to  confirm  that  our  modifications  are  correct.  Setting  pa  =  0  and 
K  =  d,  and  noting  that  p^  is  negligible  for  the  time  covered  by  the 
simulation,  we  simulate  structure  formation  in  an  Einstein  de-Sitter 
(EdS)  Universe.  We  can  use  this  freedom  to  define  three  different 
sets  of  simulation  initial  conditions. 
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Figure  1.  The  star  formation  rate  efficiency  (that  is,  star  formation  rate 
divided  by  the  total  baryon  mass  in  the  simulation  box),  for  three  simulations 
with  A  =  0  but  different  choices  for  ‘today’  [at  which  a(to)  =  1].  While 
there  is  scatter  between  the  different  simulations,  they  show  an  overall  star 
formation  history  that  is  consistent.  The  scatter  is  comparable  in  magnitude 
to  that  caused  by  using  a  different  seed  for  the  random  number  generator 
associated  with  sub-grid  physics. 


(A)  The  initial  time  of  the  simulation  has  the  same  scale 
factor  as  the  corresponding  Planck  cosmology  simulation, 
ai  j^  =  1/(1  Zi^  a)  =  1/128.  We  solve  for  the  proper  initial  time 

in  the  Planck  cosmology,  and  then  require  that  andsCO  is  normalized 
so  that  apianck(finit)  =  ^^RdsCfinit)-  This  requires  that  we  set  Hubble 
parameter  io  Ha  =  0.375. 

(B)  We  alter  the  initial  redshift  of  the  simulation  so  that  ‘today’ 
(z  =  0,  a  =  1)  is  at  ^0  =  13.8  Gyr.  This  requires  that  we  set  the  initial 
redshift  of  the  simulation  to  Zi,B  =  108  and  the  Hubble  parameter 
to  Hb  =  0.4716. 

(C)  The  Hubble  parameter  h  of  the  simulation  has  the 
same  value  as  the  corresponding  Planck  cosmology  simulation, 
^Planck  =  0.6777.  Having  found  the  time  in  the  EdS  universe  when 

=  0.6777,  we  normalize  the  scale  factor  so  that  asds  =  1  at  that 
time.  This  requires  that  we  set  the  initial  redshift  to  z/,  c  =  85.4. 

The  simulations  A,  B,  and  C  are  trying  to  solve  the  same  physical 
problem,  and  should  produce  the  same  properties  of  the  universe  as 
a  function  of  proper  time.  If  we  have  not  correctly  accounted  for 
factors  of  h  (Croton  2013)  or  confused  comoving/physical  quanti¬ 
ties  in  our  calculations,  then  these  two  simulations  should  diverge. 
Inevitably,  there  will  be  numerical  differences:  because  the  ‘time’ 
variable  of  the  simulation  is  actually  log  a,  the  time  stepping  is  not 
identical. 

Fig.  1  shows  the  star  formation  rate  efficiency  (that  is,  SFR 
divided  by  the  total  baryon  mass  in  the  simulation  box),  for  three 
simulations  (A,  B,  and  C)  with  A  =  0.  While  there  is  scatter  between 
the  different  simulations,  they  show  an  overall  star  formation  history 
that  is  consistent.  We  have  also  run  simulations  that  alter  the  seed 
for  the  random  number  generator.  The  scatter  that  this  produces  for 
a  single  set  of  parameters  is  similar  in  magnitude  to  the  differences 
between  the  simulations  A,  B,  and  C.  We  conclude  that  the  code  is 
functioning  as  expected. 

In  a  companion  paper  (Salcido  et  al.  2018),  we  consider  a  more 
detailed  comparison  between  the  EdS  cosmology  and  our  universe, 
to  quantify  the  effect  of  the  cosmological  constant  on  galaxy  for¬ 
mation  in  our  universe. 


2.4  Simulation  suite 

The  EAGLE  reference  simulations  used  cosmological  parameters 
measured  by  the  Planck  Collaboration  XVI  (2014).  We  run  seven 
EAGLE  simulations  that  modify  the  cosmological  constant,  while 
keeping  the  same  baryon  mass  per  photon  (^b).  CDM  mass  per 
photon  (^c).  and  spatial  curvature  (/c  =  0)  unchanged.  We  also  use 
the  same  physical  sub-grid  parameters  as  the  reference  model.  The 
values  of  the  cosmological  and  numerical  parameters  used  for  the 
simulations  are  listed  in  Table  2. 

As  noted  in  Section  2. 1,  we  can  solve  the  Friedmann  equations  for 
a(t)  with  an  arbitrary  normalization,  and  then  rescale  appropriately. 
For  our  cosmological  simulations,  we  choose  the  initial  scale  factor 
(or  equivalently,  redshift  Zinitiai)  to  be  the  same  for  all  values  of 
A.  In  our  Universe,  Zinitiai  =  127  corresponds  to  a  proper  time  of 
=  1 1.5  Myr.  Thus,  for  a  given  value  of  A  for  which  we  have  the 
scale  factor  a(t)  with  any  arbitrary  normalization,  we  rescale  so  that 

^(finitial)  —  1/(1  3”  -^initial)’ 

In  fact,  we  can  solve  for  the  new  cosmological  parameters  (7/q, 
^^)  in  terms  of  their  values  in  our  Universe  (Hq,  ^m) 
analytically  in  this  case.  We  require  the  expansion  of  the  universe  to 
be  the  same  at  early  times,  which  implies  that  T/g  is  equal  for  all 
universes.  In  addition,  we  increase  the  physical  energy  density  of 
dark  energy  by  a  factor/:  Anew  =/Ao,  which  implies  that  Hq  = 
-  Combining  these  equations  gives 

+  .  (10) 

^ V  +  J  ^^A 

Using  these  equations  gives  the  cosmological  parameters  in  Table  2, 
as  a  function  of  A. 

We  are  interested  in  star  formation  across  cosmic  time,  and 
so  we  want  to  run  the  simulation  as  far  as  possible  into  the  fu¬ 
ture.  This  becomes  increasingly  difficult  as  the  universe  transitions 
into  its  era  of  accelerating  expansion.  The  internal-time  variable 
in  the  code  is  log(a),  and  when  a  begins  to  increase  exponen¬ 
tially  in  cosmic  time  (/,  it  takes  more  and  more  internal-time 
steps  to  cover  the  same  amount  of  cosmic  time.  Furthermore,  be¬ 
cause  the  internal  spatial  variable  is  comoving  distance,  objects 
that  have  a  constant  proper  size  are  shrinking  in  code  units.  In 
our  experience,  in  the  accelerating  era,  the  densest  particles  in  the 
simulations  are  assigned  very  short  internal-time  steps.  The  sim¬ 
ulation  slows  to  a  crawl,  spending  inordinate  amounts  of  CPU 
time  on  a  small  number  of  particles  at  the  centres  of  isolated 
galaxies. 

In  future  work,  we  will  look  for  ways  to  overcome  these  prob¬ 
lems.  Here,  we  have  been  able  to  run  the  simulation  far  enough 
into  the  future  that,  particularly  for  large  values  of  the  cosmo¬ 
logical  constant,  quantities  such  as  the  collapse  fraction  and  the 
fraction  of  baryons  in  stars  have  approached  constant  values.  The 
endpoint  of  the  Ag  x  30,  Ag  x  100,  and  Ag  x  300  simulations 
can  be  seen  in  the  figures  in  following  sections.  We  have  captured 
the  initial  burst  of  galaxy  and  star  formation  in  these  universes, 
and  the  accelerating  expansion  of  space  makes  any  future  accre¬ 
tion  negligible.  Each  galaxy  becomes  a  separate  ‘island  universe’. 
Nevertheless,  the  far-future  ()$>20  Gyr)  fate  of  baryons  in  haloes  is 
not  captured  by  our  simulations.  Very  slow  processes  that  are  dif¬ 
ficult  to  capture  in  any  simulation  (let  alone  one  in  a  cosmological 
volume)  become  relevant:  gas  cooling  on  very  long  time-scales,  a 
trickle  of  star  formation,  rare  supernovae  in  low-density  environ¬ 
ments,  accretion  of  diffuse  gas  on  to  stellar  remnants,  and  BHs. 
These  processes  could  be  relevant  to  our  models  of  observer 
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Table  2.  Cosmological  and  numerical  parameters  for  our  simulations:  Box-size  (‘comoving’,  i.e.  the  size  of  the  box 
today  in  the  Reference  Aq  simulation),  number  of  particles,  and  cosmic  parameters  {h,  Q\),  ^^a).  (Larger  50Mpc 

boxes  were  run  and  are  analysed  in  more  detail  in  Salcido  et  al.  (2018).  For  our  purposes,  their  results  were  consistent 
with  the  25  Mpc  simulations  we  use  here.)  Note  that  these  numbers  use  the  convention  for  ‘today’  dehned  in  Section  2.4: 
the  proper  age  of  the  universe  when  ^initial  =  1/(1  +  ^initial)  is  the  same  for  all  models.  The  parameter  (7g  is  the  rms 
linear  fluctuation  in  the  mass  distribution  on  scales  of  8  Mpc,  calculated  using  camb  (Lewis  et  al.  2000).  Note  that 
this  parameter  varies  between  cosmologies  due  to  differences  in  the  growth  of  matter  fluctuations,  and  differences  in 
the  averaging  scale  due  to  the  different  values  of  h.  As  noted  in  Section  2.4,  we  do  not  run  camb  again  to  generate 
new  initial  conditions  for  each  of  our  simulations;  we  use  the  same  initial  snapshot  (particle  positions  and  velocities) 
for  each  simulation.  For  all  simulations,  the  initial  baryonic  and  dark  matter  particle  mass,  ‘comoving’  and  Plummer- 
equivalent  gravitational  softening,  and  initial  redshift  are  as  follows:  mgas  =  L81  x  10'’  Mq,  moM  =  9.70  X  lO'’  Mq, 
6com  =  2.66  kpc,  Cprop  =  0.70  kpc,  Zinitial  =  127.  Not  listed  are  the  three  simulations  used  for  the  convergence  test  (Fig.  1), 
which  use  smaller  boxes:  L=  12.5  cMpc,  N  =2  x  188^. 


Sim.  Name 

L 

[cMpc] 

N 

h 

^A 

0^8 

EdS_25  (A  =  0) 

25 

2 

X  376^ 

0.3755 

1.0 

0.1572 

0 

0.6826 

Ref_25  (Ao) 

25 

2 

X  376^ 

0.6777 

0.307 

0.048  25 

0.693 

0.8288 

Aq  X  3 

25 

2 

X  376^ 

1.047 

0.1287 

0.0202 

0.8713 

0.8913 

Ao  X  10 

25 

2 

X  376^ 

1.823 

0.0424 

0.006  67 

0.9576 

0.8955 

Ao  X  30 

25 

2 

X  376^ 

3.113 

0.01455 

0.002  29 

0.985  45 

0.8434 

Ao  X  100 

25 

2 

X  376^ 

5.654 

0.00441 

6.93  X  10"^ 

0.995  59 

0.7476 

Ao  X  300 

25 

2 

X  376^ 

9.779 

0.001  47 

2.32x  lO"'* 

0.998  53 

0.6446 

creation  over  all  of  cosmic  time;  we  will  return  to  these  issues  in 
Section  5.3. 


3  CHANGING  THE  COSMOLOGICAL 
CONSTANT:  GLOBAL  PROPERTIES 

We  vary  the  cosmological  constant  between  zero  and  several  hun¬ 
dred  times  larger  than  the  value  in  our  Universe.  We  do  not  consider 
negative  values  of  the  cosmological  constant  here,  as  it  would  re¬ 
quire  significant  changes  to  the  time  stepping  in  the  code  to  handle 
the  transition  from  expansion  to  contraction. 

Fig.  2  shows  the  deceleration  parameter  q  =  —a/(aH^)  and  the 
linear  growth  factor  D(t)  as  a  function  of  cosmic  time,  for  different 
values  of  the  cosmological  constant.  As  the  cosmological  constant 
increases,  the  time  at  which  the  expansion  of  the  universe  begins  to 
accelerate  (q  <  0)  moves  to  earlier  times  as  4ccei  1  /-y/GpA-  Once 
accelerated  expansion  begins,  the  formation  of  structure  freezes 
and  accretion  stops.  We  can  see  this  in  linear  perturbation  theory, 
where  all  modes  grow  in  proportion  to  the  growth  factor  D(t);  we 
normalize  D(t)  so  the  curves  are  equal  at  early  times,  and  D(to)  =  1 
in  our  Universe  today.  We  see  that  once  the  expansion  of  the  universe 
begins  to  accelerate,  the  growth  factor  approaches  a  constant,  and 
structures  ceases  to  grow. 

In  this  section,  we  will  characterize  the  details  of  structure  for¬ 
mation  in  these  universes.  Ordinarily,  one  describes  these  proper¬ 
ties  using  comoving  quantities,  such  as  the  comoving  halo  number 
density  and  comoving  star  formation  rate  density.  One  immedi¬ 
ate  problem  is  that  the  term  ‘comoving’  is  meaningless  when  dif¬ 
ferent  universes  are  being  compared.  There  is  no  ‘today’  that  is 
common  to  all  models,  relative  to  which  we  can  define  comov¬ 
ing  volumes,  densities,  and  the  like.  There  is  nothing  special, 
cosmically  speaking,  about  13.8  Gyr  or  2.725  K.  We  can  arbi¬ 
trarily  change  the  comoving  density  of  star  formation,  for  exam¬ 
ple,  by  choosing  a  different  cosmic  time  in  a  given  universe  to 
be  ‘today’,  which  makes  the  comparison  of  comoving  densities 
meaningless. 


Figure  2.  The  deceleration  parameter  q  =  —a/{aH^)  (top)  and  the  linear 
growth  factor  D(t)  (bottom)  as  a  function  of  cosmic  time,  for  different  values 
of  the  cosmological  constant.  Note  that  q=  1  /2  at  all  times  for  the  A  =  0 
cosmology.  As  the  cosmological  constant  increases,  the  time  at  which  the 
expansion  of  the  universe  begins  to  accelerate  (q  <  0)  moves  to  earlier  times 
as  7accei  ~  l/\/UpA-  Once  accelerated  expansion  begins,  the  formation  of 
structure  freezes  and  accretion  stops,  and  D(t)  approaches  a  constant. 
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Figure  3.  Left:  the  fraction  of  mass  in  each  simulation  that  is  part  of  a  resolved  halo:  F(  >  mminlO>  where  mmin  =  3.1  x  10*  Mq.  This  minimum  mass  is  a 
consequence  of  the  numerical  resolution  of  the  simulations,  but  is  consistent  across  all  of  them  and  approximately  excludes  haloes  that  are  too  small  to  form 
stars.  The  result  is  a  measure  of  the  fraction  of  mass  in  the  universe  that  resides  in  potentially  star-forming  haloes.  Right:  the  specific  total  halo  accretion  rate, 
i.e.  the  time  derivative  of  the  left  hand  curve.  The  rate  peaks  dXt  =  0.8  Gyr  in  our  Universe  (A  =  Aq).  Even  for  a  universe  with  a  cosmological  constant  10 
times  larger  than  ours  (Aq  x  10),  there  is  a  minimal  difference  in  total  halo  mass  fraction  even  after  20  Gyr,  well  into  the  accelerating  phase  of  the  universe’s 
expansion. 


To  overcome  this,  we  will  calculate  quantities  relative  to  the  phys¬ 
ical  mass  (total  or  baryonic)  in  the  simulation  box.^  This  provides  a 
meaningful  comparison  between  the  simulated  universes,  and  like 
comoving  densities  it  does  not  automatically  scale  with  expansion 
of  the  universe.  We  can  ask,  for  example,  what  fraction  of  the  total 
baryonic  mass  in  the  universe  is  in  the  form  of  stars  as  a  function  of 
cosmic  time?  What  fraction  has  been  converted  into  metals? 


3.1  Mass  accretion 

Formally,  in  a  CDM  universe,  every  particle  is  in  a  dark  matter  halo 
of  some  mass.  That  is,  the  collapse  fraction  of  the  universe  is  always 
unity;  from  Press  &  Schechter  (1974)  theory, 

F(>M|0  =  erfc  (H) 

VV2a(M,  t)J 

where  F{>M\t)  is  the  fraction  of  matter  at  cosmic  time  t  that  is  part 
of  a  collapsed  halo  of  mass  greater  than  M,  <5crit(0  is  the  critical 
linear  overdensity  of  a  collapsed  object,  and  cr(M,  i)  is  the  standard 
deviation  of  the  cosmic  matter  field  when  smoothed  on  a  scale  that 
encloses  mass  M.  The  matter  variance  <7(M,  0  ^  oo  as  M  ^  0, 
thus  F(>0|0  =  1  at  all  times. 

In  the  simulation,  however,  there  is  a  minimum  dark  matter  halo 
mass  that  can  be  resolved  by  the  particles.  Given  that  each  dark 
matter  particle  has  mass  moM  =  9.7  x  10^  Mq  and  we  require 
32  particles  to  identify  a  halo,  we  can  resolve  haloes  with  mass 
greater  than  rnmin  =  3.1  x  10^  Mq.  Summing  the  total  mass  in 
these  haloes,  then,  gives  the  collapse  fraction  for  resolved  haloes: 
F(>mminl0-  This  approximately  excludes  haloes  that  are  too  small 
to  form  stars,  so  gives  us  the  fraction  of  mass  in  the  universe  that 
resides  in  potentially  star-forming  haloes;  the  remainder  can  be 
considered  as  the  intergalactic  medium. 

^  We  could  define  comoving  densities  relative  to  the  initial  physical  volume 
of  the  simulation  box,  which  is  the  same  in  all  models.  But  while  this  allows 
a  meaningful  comparison,  the  initial  cosmic  time  is  still  arbitrary.  Choosing 
an  earlier  time  would  increase  all  the  ‘comoving’  densities,  which  makes 
their  value  in  a  given  universe  difficult  to  interpret.  Calculating  specific  (per 
unit  mass)  quantities  overcomes  this  problem. 


Fig.  3  shows  (left)  the  fraction  of  the  total  mass  in  the  universe 
that  has  collapsed  into  resolved  haloes,  and  {right)  the  specific  total 
halo  accretion  rate,  i.e.  the  time  derivative  of  the  left-hand  curve. 
In  this  figure  and  those  following,  the  time  derivative  is  calculated 
after  smoothing  the  accretion  fraction.  Even  for  a  universe  with  a 
cosmological  constant  10  times  larger  than  ours  (Aq  x  10),  there  is 
minimal  difference  in  total  halo  mass  fraction  even  after  20  Gyr,  well 
into  the  accelerating  phase  of  the  universe’s  expansion.  The  initial 
peak  in  the  accretion  rate  at  ?  =  0.8  Gyr  remains  largely  unchanged 
even  in  a  universe  with  a  cosmological  constant  30-100  times  larger 
than  ours.  In  a  universe  with  A  =  Aq  x  100,  a  fifth  of  the  mass  in 
the  universe  accretes  into  haloes  with  m  >  rnmin  =  3.1  x  10^  Mq. 


3.2  Baryon  flow 

Baryons  are  subject  to  physical  forces  other  than  gravity:  the 
smoothing  effect  of  gas  pressure,  cooling  and  heating  from  radi¬ 
ation,  star  formation,  supernovae-driven  galactic  winds,  BH  feed¬ 
back,  and  more.  Fig.  4  shows,  left,  the  fraction  of  the  baryonic  mass 
(in  the  form  of  stars  and  gas)  in  the  simulation  that  is  inside  dark 
matter  haloes  with  m  >  rnmin  =  3.1  x  10^  Mq  as  a  function  of 
cosmic  time,  and,  right,  the  specific  rate  of  baryon  accretion  (i.e. 
per  unit  total  baryon  mass). 

We  see  the  same  peak  in  the  accretion  rate  at  t  =  0.8  Gyr,  and 
when  there  is  zero  cosmological  constant,  the  baryon  accretion  rate 
increases  in  a  similar  way  to  the  total  accretion  rate  (Fig.  3).  As  the 
cosmological  constant  increases,  it  has  a  much  larger  effect  on  the 
baryons  than  the  dark  matter.  In  fact,  for  A  =  Aq  x  10,  the  rate  of 
baryon  accretion  becomes  negative,  as  baryons  are  -  on  average  - 
being  ejected  from  galaxies. 

We  can  understand  this  effect  as  follows.  We  can  write  the  accel¬ 
eration  (Og)  of  a  test  mass  at  distance  r  from  a  large  mass  M  under  the 
Newtonian  gravitational  force  with  a  cosmological  constant  term. 


M  Ac^ 

a,  =  -G-  +  — r. 


(12) 


If  we  consider  a  large  collapsed  mass,  then  the  distance  (do)  at 
which  the  force  on  a  test  mass  is  balanced  between  attraction  to  the 
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Figure  4.  Left  :  the  fraction  of  the  baryonic  mass  in  the  simulation  that  is  inside  dark  matter  haloes  with  m  >  nimin  =  3.1  x  10*  Mq  as  a  function  of  cosmic 
time,  and  Right:  the  specific  rate  of  baryon  accretion.  The  rate  peaks  at  t  =  0.6  Gyr  in  our  Universe  (A  =  Aq). 


Figure  5.  Left:  the  fraction  of  cosmic  baryons  by  mass  that  are  in  the  form  of  stars  as  a  function  of  cosmic  time.  Right:  the  star  formation  rate  efficiency, 
which  takes  into  account  both  star  birth  only.  The  rate  peaks  at  t  =  3.5  Gyr  in  our  Universe  (A  =  Aq). 


central  mass  and  repulsion  by  the  cosmological  constant  is  found 
by  setting  Ug  =  0, 

/  M  \  A  \ 

U)  ’ 

or  equivalently  in  terms  of  the  ratio  Pa/Paq-  our  universe,  this  is 
^4  times  larger  than  the  virial  radius  of  the  halo  (which  also  scales 
as  the  1/3  power  of  mass).  In  universes  in  which  the  cosmological 
constant  is  larger,  these  distances  are  comparable. 

As  seen  in  Fig.  4,  this  does  not  dramatically  affect  the  growth  of 
the  dark  matter  halo.  But  baryonic  matter  ejected  from  galaxies  in 
galactic  winds  or  outflows,  if  it  reaches  the  outer  parts  of  the  halo, 
is  liable  to  be  lost.  Rather  than  raining  back  down  on  the  galaxy 
after  a  delay  of  '^l  Gyr  (Oppenheimer  &  Dave  2008;  Oppenheimer 
et  al.  2010;  van  de  Voort  2017),  this  material  is  lost,  drawn  away 
into  the  expansion  of  the  universe  by  the  repulsive  effect  of  the 
cosmological  constant  (Barnes  et  al.  2006). 

As  we  will  see  in  the  next  subsection,  in  universes  for  which 
(A  Ao  X  10),  the  initial  burst  of  star  formation  in  the  universe 
occurs  when  the  universe  has  begun  to  expand  exponentially.  This 
rapid  star  formation,  combined  with  BH  feedback,  launches  out¬ 


flows  that  are  carried  away  by  the  accelerating  expansion.  This 
effect  overwhelms  accretion  by  gravitational  attraction,  causing  the 
net  accretion  rate  to  become  negative.  The  result  is  that  there  is 
not  a  simple,  linear  relationship  between  dark  matter  halo  growth 
and  baryon  accretion  that  holds  for  all  values  of  the  cosmological 
constant. 

3.3  Star  formation 

Some  of  the  baryons  that  accrete  into  haloes  will  form  stars.  Fig.  5 
shows  {left)  the  fraction  of  cosmic  baryons  that  are  in  the  form  of 
stars  as  a  function  of  cosmic  time,  and  {right)  the  star  formation  rate 
efficiency,  which  takes  into  account  star  birth  only.  Note  that,  as  it 
is  commonly  used  in  the  galaxy  formation  literature,  ‘speciflc  star 
formation’  refers  to  the  star  formation  rate  of  a  galaxy  divided  by 
its  stellar  mass.  To  avoid  confusion,  we  will  call  the  star  formation 
mass  (rate)  per  unit  total  baryon  mass  the  star  formation  (rate) 
efficiency. 

The  star  formation  rate  efficiency  peaks  at  t  ~  3.5  Gyr.  This  is 
a  delayed  consequence  of  the  peak  in  the  mass  accretion  rate  at 
t  =  0.8  Gyr.  As  the  cosmological  constant  increases,  the  haloes 
are  starved  both  by  the  cessation  of  fresh  accretion  from  the 
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Figure  6.  Left:  the  fraction  of  cosmic  baryons  that  are  in  the  form  of  metals  in  collapsed  haloes.  Right:  the  halo  metal  production  rate.  The  rate  peaks  at 
t  =  3.2  Gyr  in  our  Universe  (A  =  Aq),  and  the  peak  is  steadily  diminished  as  the  cosmological  constant  increases. 


intergalactic  medium  and  by  the  lack  of  recycling  of  outflowing 
gas,  noted  above.  The  result  is  a  signiflcant  curtailing  of  the  star 
formation  rate  efficiency.  While  the  A  =  0  universe  has  turned 
~4  per  cent  of  its  baryons  into  stars  by  ?  =  20  Gyr,  for  Aq  x  100, 
this  fraction  is  essentially  constant  after  10  Gyr  at  0.5  per  cent.  This 
factor  of  8  decrease  contrasts  with  the  factor  of  2.4  decrease  in  the 
total  mass  accretion. 

3.4  Cosmic  metal  production 

Planets  and  their  occupants  are  formed  from  the  products  of  stellar 
nucleosynthesis.  The  eagle  code,  in  addition  to  primordial  H  and 
He,  follows  nine  metals:  C,  N,  O,  Ne,  Mg,  Si,  S,  Ca,  and  Fe.  Fig.  6 
shows  (left)  the  fraction  of  cosmic  baryons  that  are  in  the  form  of 
metals  in  collapsed  haloes,  and  {right)  the  halo  metal  production 
rate.  Note  that  this  includes  metals  in  all  phases:  inside  stars,  in  dense 
star-forming  clouds,  and  in  the  hot,  non- star  forming  interstellar  gas. 
The  halo  metal  production  rate  reflects  the  balance  between  metal 
production  in  stars,  recycling  back  into  the  interstellar  medium  by 
winds  and  supernovae,  re-incorporation  into  later  generation  stars, 
ejection  from  haloes  in  galactic  winds,  and  reaccretion  into  haloes. 

As  star  formation  peaks  (Fig.  5),  metals  are  being  produced  in 
stars  and  returned  to  the  IGM  in  supernovae  and  planetary  nebulae. 
This  feedback  also  loads  metals  into  the  galactic  winds  that  drive 
baryons  out  of  haloes  (Fig.  4).  As  with  the  baryon  accretion  rate,  the 
net  accretion  rate  becomes  negative  for  certain  values  of  A  as  metals 
are  ejected  in  winds  at  a  higher  rate  than  they  are  produced  and  reac¬ 
creted.  Our  universe  turns  approximately  a  fraction  1.2  x  10  ^  by 
mass  of  its  baryons  into  halo  metals  by  20  Gyr,  while  for  Aq  x  100 
the  fraction  asymptotes  by  10  Gyr  to  1.5  x  10  This  factor  of  10 
difference  contrasts  with  the  factor  of  2.5  difference  with  regards  to 
the  total  fraction  of  mass  in  haloes. 

4  ACCRETION  AND  STAR  FORMATION  IN 
INDIVIDUAL  HALOES 

In  this  section,  we  will  consider  the  evolution  of  a  comoving  region 
of  space  that,  in  our  Universe,  evolves  into  a  Milky  Way-mass  halo 
(2  X  10^^  Mq)  by  the  present  day.  Fig.  7  shows  the  projected  gas 
density  in  a  comoving  region  around  the  halo  equal  to  4  Mpc  today 
in  our  universe;  the  cosmic  time  and  proper  size  of  the  region  are 
shown  above  each  panel.  The  /^^-hand  panels  show  an  EdS  universe 
(A  =  0);  the  r/g/tUhand  panels  show  a  Aq  x  30  universe. 


The  top  two  panels  show  this  region  of  the  universe  at  cosmic 
time  t  =  0.757  Gyr,  while  the  Aq  x  30  is  still  in  its  early  decelerating 
phase.  The  proper  sizes  of  the  boxes  are  within  1  percent  of  each 
other,  and  the  distributions  of  matter  are  very  similar.  We  see  the 
usual  picture  of  small  haloes  collapsing  and  hierarchically  merging 
into  larger  haloes. 

The  middle  two  panels  show  this  region  of  the  universe  at  cosmic 
time  t  =  6.5  Gyr.  The  Aq  x  30  is  undergoing  accelerating  expansion, 
so  the  proper  size  of  the  region  is  2.3  times  larger  than  in  the  EdS 
universe,  and  the  linear  growth  factor  is  33  per  cent  smaller.  The 
large  central  halo  in  the  EdS  simulation  has  drawn  in  a  more  matter 
from  its  surroundings,  and  is  still  being  drawn  towards  a  second 
halo  at  the  bottom  of  the  panel. 

The  bottom  two  panels  show  this  region  of  the  universe  at  cosmic 
time  t  =  12.5  Gyr.  The  accelerating  expansion  of  the  Aq  x  30 
means  that  the  proper  size  of  the  comoving  region  is  10  times  larger 
than  in  the  EdS  universe.  The  typical  Newtonian  force  between 
two  masses  in  the  region  is  thus  100  times  smaller,  and  the  linear 
growth  factor  is  2.3  times  smaller.  The  difference  in  the  distribution 
of  matter  is  quite  dramatic:  in  the  Aq  x  30  universe,  there  has  been 
little  evolution  of  the  structure  of  the  universe  since  t  =  6.5  Gyr. 
The  matter  in  the  vertical  filaments  has  not  fallen  into  the  large 
halo,  starving  the  galaxy  of  gas.  In  the  EdS  universe,  the  halo  has 
been  drawn  closer  to  the  second  halo  at  the  bottom  of  the  panel;  the 
filament  of  matter  between  them  has  largely  fallen  into  one  of  the 
haloes. 

To  highlight  the  difference  between  the  final  states  of  the  galaxies 
at  the  centre  of  the  halo,  Eig.  8  shows  a  region  of  constant  proper 
size  (0.5  Mpc)  around  the  central  galaxy  in  the  regions  shown  in 
Eig.  8.  The  colour  scaling  in  all  four  panels  is  held  constant. 

The  top  two  panels  show  this  region  of  the  universe  at  cosmic 
time  t  =  6.5  Gyr.  Both  show  a  galaxy  in  formation,  being  fed  by 
streams  of  gas.  But  already  we  can  see  that  the  EdS  galaxy  {left)  is 
larger,  and  is  surrounded  by  a  much  higher  density  circumgalactic 
medium.  In  the  Aq  x  30  universe  {right),  the  free  fall  time  from 
the  edge  of  the  isolated  region  around  the  galaxy  is  a  few  Gyr,  and 
so  the  halo  accretes  as  much  material  as  is  available  on  this  time- 
scale.  Accordingly,  the  total  mass  of  the  halo  only  grows  by  only 
~  1  per  cent  in  the  6  Gyr  between  the  two  snapshots  shown  in  Eig.  8, 
to  a  final  mass  of  8  x  10^^  Mq.  In  this  isolation,  the  gas  collapses 
into  one  monolithic  disc.  In  the  same  time  in  the  EdS  simulation, 
the  halo  has  doubled  in  mass  to  4  x  10^^  Mq  at  ^  =  12.5  Gyr,  and 
is  still  growing. 
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EdS:  t  =  0.757  Gyr,  box  =  0.497  pMpc  Lambda  x  30:  t  =  0.757  Gyr,  box  =  0.504  pMpc 


EdS:  t  =  12.465  Gyr,  box  =  3.208  pMpc  Lambda  x  30:  t  =  12.465  Gyr,  box  =  3 1 .790  pMpc 


Figure  7.  The  evolution  of  the  projected  gas  density  of  a  comoving  region  of  space  that,  in  our  universe,  evolves  into  a  Milky  Way-mass  halo  by  the  present 
day.  The  comoving  size  is  4Mpc  in  our  Universe.  The  proper  time  and  proper  size  of  the  region  are  shown  above  each  panel.  Left:  an  EdS  universe  (A  =  0); 
Right:  a  Ao  X  30  universe.  The  colour  scaling  on  each  row  is  held  constant. 
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EdS:  t  =  6.468  Gyr,  box  =  0.5  pMpc  Lambda  x  30:  t  =  6.468  Gyr,  box  =  0.5  pMpc 


Figure  8.  The  evolution  of  the  projected  proper  gas  density  in  a  region  of  fixed  proper  size  (0.5  Mpc)  around  the  central  galaxies  in  the  panels  in  Fig.  8.  The 
proper  time  is  shown  above  each  panel.  Left:  an  EdS  universe  (A  =  0);  Right:  a  Aq  x  30  universe. 


5  IMPLICATIONS  FOR  MULTIVERSE  MODELS 
5.1  The  measure  of  the  multi  verse 

We  can  use  our  calculations  to  make  predictions  from  multiverse 
models.  Given  a  model  that  predicts  an  ensemble  of  universes  with 
a  distribution  of  values  for  the  cosmological  constant,  we  can  ask 
what  fraction  of  observers  will  inhabit  a  universe  with  a  particular 
value  of  A. 

If  the  model  in  question  predicts  definite  ensemble  of  universes, 
inhabited  by  a  finite  number  of  observers,  then  this  calculation  is 
straightforward.  Scientific  theories  are  tested  by  predicting  obser¬ 
vations,  and  so  all  observers  are  treated  as  of  equal  importance  for 
the  purposes  of  calculating  the  likelihood.^  We  thus  use  a  counting 


^We  will  ignore  the  complication  of  asking:  what  exactly  counts  as  an 
observer?  We  cannot  predict  the  occurrence  of  observers  in  sufficient  detail 
to  make  any  difference.  That  is,  we  might  wonder  whether  any  complex  life 
form  counts  as  an  observer  (an  ant?),  or  whether  we  need  to  see  evidence  of 
communication  (a  dolphin?),  or  active  observation  of  the  universe  at  large 
(an  astronomer?).  Our  model  does  not  contain  anything  as  detailed  as  ants, 


metric  to  calculate  the  likelihood, 

/  K  \  Ti^  px  1  .  ^obs([A,  A  +  dA]) 

/7obs(A|M5)dA  =  - ,  (14) 

^obs 

where  M  is  the  multiverse  model,  B  is  any  relevant  background 
information  (which  should  not  give  away  any  clues  about  the  prop¬ 
erties  of  the  actual  universe),  7tobs([A,  A  +  dA])  is  the  number  of 
observers  (or  observer-moments)  that  exist  in  a  universe  with  cos¬ 
mological  constant  in  the  range  [A,  A  +  dA],  and  TZobs  is  the  total 
number  of  observers  in  the  multiverse. 

To  evaluate  these  quantities,  we  calculate  (at  least  approximately) 
the  rate  at  which  observers  are  produced  per  unit  time  per  unit  co¬ 
moving  volume,  for  a  given  set  of  cosmic  and  fundamental  parame¬ 
ters:  d^TZobs/dMV.  So  long  as  the  universe  has  finite  age,  or  if  the  rate 
at  which  observers  are  produced  approaches  zero  quickly  enough 
into  the  future,  then  the  integral  over  cosmic  time  of  this  rate  will 


dolphins  or  astronomers,  so  we  are  unable  to  make  such  a  fine  distinction 
anyway.  In  any  case,  such  a  distinction  is  unlikely  to  bias  our  calculation 
towards  any  particular  value  of  the  cosmological  constant. 
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be  finite.  Then,  the  likelihood  of  the  cosmological  constant  is 


Pobs(A|MB)dA 


j;-""  V(f;A)  pvim^dtdA 
V(t-  A)  pvim^dtdA  ’ 


(15) 


where  ?max  is  the  maximum  age  of  the  universe  (possibly  infinite), 
V(t;  A)  is  the  total  comoving  volume  of  the  universe,  pv(A\t)dA  is 
the  fraction  of  the  universe  by  comoving  volume  at  time  t  in  which 
the  value  of  the  cosmological  constant  is  in  the  range  [A,  A  +  dA]. 
The  comoving  volume  depends  on  the  arbitrary  normalization  of 
a(t),  but  this  cancels  in  the  equation  above. 

However,  most  proposed  multiverses  are  not  finite.  In  eternally 
inflating  universes,  for  example,  it  is  argued  that  not  only  does 
the  multiverse  consist  of  an  infinite  number  of  universes,  but  each 
universe  is  infinitely  large  (Vilenkin  &  Winitzki  1997;  Garriga  & 
Vilenkin  2001;  Knobe,  Olum  &  Vilenkin  2003;  Freivogel  et  al. 
2006;  Guth  2007;  Ellis  &  Stoeger  2009).  Thus,  the  number  of  uni¬ 
verses  with  a  given  value  of  A,  times  the  average  number  of  ob¬ 
servers  in  those  universes,  divided  by  the  total  number  of  observers 
in  the  multiverse,  is  oo  x  oo/oo. 

These  infinities  need  to  be  managed  with  a  measure;  see,  among 
others,  Vilenkin  (1995),  Garriga  et  al.  (2006),  Aguirre,  Gratton 
&  Johnson  (2007),  Vilenkin  (2007a),  Vilenkin  (2007b),  Gibbons 
&  Turok  (2008),  Page  (2008),  Bousso,  Freivogel  &  Yang  (2009), 
de  Simone  et  al.  (2010),  Freivogel  (2011),  Bousso  &  Susskind 
(2012),  Garriga  &  Vilenkin  (2013),  and  Page  (2017).  Simplistically, 
this  measure  can  be  thought  of  in  two  ways.  First,  a  multiverse 
model  could  motivate  confining  our  attention  to  a  finite  region  of 
the  universe  with  volume  V(t;  A)  (as  a  function  of  time  and  A).  Then, 
we  can  use  the  finite  calculation  for  the  likelihood  (equation  15). 
Secondly,  the  measure  could  specify  the  fraction  of  the  volume  of 
the  universe  in  which  cosmic  parameters  are  in  a  given  range,  even 
though  the  total  volume  of  the  universe  is  infinite.  This  is  used  to 
weight  the  integral,  effectively  ‘cancelling’  the  infinite  quantity  V(t; 
A)  from  the  numerator  and  denominator  of  equation  (15),  which 
gives 


PoUA\MB)dA 


/o"'"Bv(A|o|t^dfdA 
JJf- py(A\t)^dtdA' 


(16) 


Here,  rather  than  focus  on  a  specific  multiverse  model,  we  will 
consider  three  measures.  Following  Weinberg  (1987),  Efstathiou 
(1995),  Peacock  (2007),  and  Bousso  &  Leichenauer  (2010),  we 
note  that  nothing  in  fundamental  physics  picks  out  a  value  of  the 
cosmological  constant  as  privileged,  including  the  value  zero.  This, 
in  particular,  rules  out  the  use  of  a  logarithmic  prior.  In  the  range 
of  A  that  we  consider,  which  is  very  small  compared  to  the  Planck 
scale,  we  approximate  the  distribution  as  flat  on  a  linear  scale.  The 
difference  between  the  measures  is  the  quantity  with  respect  to 
which  the  distribution  is  flat. 


1.  Mass-weighted:  there  is  a  uniform  probability  that  a  given 
mass  element  in  the  universe  will  inhabit  a  region  with  a  given 
value  of  the  cosmological  constant.  Note  that,  for  reasons  discussed 
in  Section  3,  specifying  that  there  is  uniform  probability  with  respect 
to  comoving  volume  is  not  sufficient,  as  there  is  no  universal  ‘today’ 
relative  to  which  we  can  define  volume.^  We  use  the  constraint  of 
constant  mass  to  define  comoving  volumes  between  universes. 


^  Put  another  way,  we  are  free  to  renormalize  a{t),  but  this  normalization 
could  depend  on  A.  This  will  not  cancel  out  in  equation  (16),  making  the 
calculated  probability  arbitrary. 


2.  Causal  patch:  this  measure  was  proposed  to  solve  the  quantum 
xeroxing  paradox  in  BHs  (Susskind,  Thorlacius  &  Uglum  1993; 
Bousso  2006;  Bousso,  Freivogel  &  Yang  2006),  treating  the  de- 
Sitter  horizon  in  a  universe  with  A  analogously  to  a  BH  horizon. 
We  ask:  what  is  the  volume  of  the  region  of  the  universe  at  time  t 
that  can  causally  affect  a  given  comoving  world  line  in  the  future 
of  tl  The  comoving  extent  of  the  region  is 

/•w 

Xpatch(0  ~  /  V7\’  (1"^) 

Jt  a{t) 

Then,  for  the  spatially  flat  universes  that  we  consider  here,  the 
volume  is  y  =  (47r/3)x^(0,  which  goes  in  equation  (15).  Note  that 
equation  (17)  depends  on  the  arbitrary  normalization  of  a{t),  but 
this  is  cancelled  out  when  it  is  multiplied  by  the  observer  creation 
rate:  d^^Uou/dtdV.  The  comoving  size  of  the  causal  patch  is  shown  in 
Fig.  9  {left,  relative  to  the  normalization  of  a{i)  from  Section  2.1). 
Also  shown  {middle)  is  the  physical  mass  contained  within  the 
causal  patch  (which  is  not  relative  to  the  normalization  of  a{i)). 

3.  Causal  diamond:  this  measure  is  based  on  the  principle  that 
space-time  regions  that  are  causally  inaccessible  should  be  disre¬ 
garded  (Bousso  2006;  Bousso  et  al.  2007).  We  consider  a  comoving 
world  line  in  a  universe,  extending  from  the  end  of  inflation  (reheat¬ 
ing  at  ^  =  ^rh)  to  the  distant  future.  What  is  the  volume  (at  time  t)  of 
the  region  of  the  universe  that  is  enclosed  by  a  photon  that  departs 
the  world  line  at  its  beginning  and  returns  at  the  end?  We  can  write, 

Xdiamond(0  —  nai^{Xpatch(0?  ^(01?  (1^) 

where  r]{t) 

As  for  the  causal  patch,  the  volume  V  =  (47r/3)x^(0  is  used  in 
equation  (15).  The  causal  diamond  is  shown  in  Fig.  9.  Also  shown 
is  the  physical  mass  contained  within  the  causal  diamond  {right). 

We  stress,  however,  that  the  measure  is  not  a  ‘degree  of  free¬ 
dom’  in  a  multiverse  model.  It  must  not  be  inferred  from  or  fit 
to  observations,  and  the  fact  that  a  particular  measure  gives  good 
agreement  with  observations  is  no  reason  to  prefer  that  measure. 
The  reason  is  that  any  value  of  A  can  be  made  practically  certain 
with  an  appropriately  jerry-rigged  measure.  If  a  model  derives  its 
prediction  from  observations,  then  its  predictions  cannot  then  be 
tested  by  those  same  observations.  A  multiverse  model  is  supposed 
to  tell  us  about  the  global  structure  of  the  universe.  There  should  not 
be  any  assumptions  that  need  to  be  added  ‘on  top’ ,  because  there  are 
no  physical  facts  left  to  specify,  at  least  on  relevant  cosmological 
scales.  The  measure  should  follow  naturally  -  in  some  sense  -  from 
the  multiverse  model. ^ 

5.2  Models  of  observers 

We  need  to  connect  the  presence  of  observers  to  local  conditions  in 
our  simulations.  This  will,  inevitably,  be  a  combination  of  approxi¬ 
mation  and  guesswork.  Note  that  any  constant  factor  in  the  observer 

^  To  put  this  another  way,  suppose  a  multiverse  model  specified  the  global 
structure  of  the  universe  in  painstaking  detail:  the  value  of  cosmic  parameters 
and  properties  at  every  place  and  time.  What  would  it  mean  to  apply  two 
different  measures  to  this  model,  to  derive  two  different  predictions?  How 
could  all  the  physical  facts  be  the  same,  and  yet  the  predictions  of  the  model 
be  different  in  the  two  cases?  What  is  the  measure  about,  if  not  the  universe? 
Is  it  just  our  own  subjective  opinion?  In  that  case,  you  can  save  yourself 
all  the  bother  of  calculating  probabilities  by  having  an  opinion  about  your 
multiverse  model  directly. 


dt 

t,h  a(0 


(19) 
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Figure  9.  Left:  comoving  causal  patch  and  causal  diamond  versus  proper  time  (Gyrs)  for  different  values  of  the  cosmological  constant  shown  in  the  legend 
of  Fig.  11.  The  decreasing  curves  are  the  causal  patch.  Increasing  (overlapping)  curves  are  the  quantity  v(t)  from  equation  (18);  the  causal  diamond  is  the 
minimum  of  these  two  curves  at  a  given  time.  The  comoving  distance  is  relative  to  the  chosen  normalization  of  a(t),  as  described  in  Section  2.1.  Also  shown 
are  the  physical  mass  inside  the  causal  patch  {middle)  and  causal  diamond  {right)  as  a  function  of  cosmic  time,  which  are  independent  of  the  normalization. 


creation  rate  will  cancel  in  equations  (15)  and  (16),  so  an  absolute 
rate  is  not  required.  We  consider  three  models  of  observers,  linked 
to  the  production  of  energy  and  chemical  elements. 


(1)  Star  formation  +  fixed  delay:  following  Bousso  &  Le- 
ichenauer  (2010),  we  consider  a  model  in  which  observers  follow 
the  formation  of  a  star  with  a  fixed  time  delay  of  5  Gyr.  We  also 
considered  a  time  delay  of  10  Gyr,  but  it  made  minimal  difference 
to  our  conclusions.  This  is  inspired  by  the  time  taken  for  intelligent 
life  to  form  on  Earth  after  the  birth  of  the  Sun. 

(2)  Star  formation  +  main-sequence  lifetime:  as  first  argued  by 
Carter  (1983,  see  also  Barrow  &  Tipler  1986),  if  the  formation  of 
life  is  extremely  improbable  -  that  is,  if  the  average  time-scale  for 
its  formation  is  much  longer  than  the  lifetimes  of  stars  -  then  it 
will  form  at  the  last  available  moment,  so  to  speak.  Most  stars  will 
host  lifeless  planets,  but  where  life  forms  it  will  do  so  at  a  time  that 
is  of  the  order  of  the  main- sequence  lifetime  of  the  star.  As  a  first 
approximation,  we  assume  that  there  is  a  constant  probability  per 
unit  time  of  life  forming  around  stars  of  all  masses.  The  observer 
creation  rate  for  each  star  population  that  forms  is  proportional  to 
the  fraction  of  stars  (by  number)  that  are  still  on  the  main  sequence 
after  time  A?, 


/ms(A0 


J  e(t^,{M)  -  At)  g(M)dM 
/  t(M)dM 


(20) 


where  ^(M)  is  the  stellar  initial  mass  function  (IMF),  tra^{M)  is  the 
main- sequence  lifetime  of  a  star  of  mass  M,  the  limits  of  the  in¬ 
tegral  are  the  minimum  and  maximum  stellar  masses,  and  0{x)  is 
the  Heaviside  step  function,  so  that  only  those  stars  whose  main- 
sequence  lifetimes  are  longer  than  the  time  since  the  population  was 
born  contribute.  We  use  the  Chabrier  (2003)  IMF,  and  a  simple  rela¬ 
tionship  between  mass  and  main- sequence  lifetime  drawn  from  the 
analytic  model  of  Adams  (2008)  normalized  to  tms  =  10  Gyr  at  solar 
mass;  this  is  broadly  consistent  with  Portinari,  Chiosi  &  Bressan 
(1998).  Of  particular  importance  are  the  maximum  and  minimum 
stellar  masses.  To  be  consistent  with  the  IMF  used  to  calibrate  the 
EAGLE  simulations,  we  choose  the  minimum  and  maximum  stellar 
masses  to  be  Mmin  =  0.1  Mq  and  M^ax  =  100  Mq.  The  resulting 
main- sequence  fraction  is  shown  in  Fig.  10. 

Folding  in  the  star  formation  (birth)  rate  density  (Pstar),  we  calculate 
the  global  observer  creation  rate.  A  stellar  population  that  formed 
at  time  At  before  the  present  time  t  provides  a  relative  contribution 
of/ms(At)  to  the  observer  creation  rate 


dtdV 


(0  oc 


/>star(^')/ms(^  -  t')&d 


(21) 


Figure  10.  The  fraction  of  stars  by  number  that  are  still  on  the  main  se¬ 
quence  of  their  evolution  after  time  At.  We  assume  a  Chabrier  (2003)  IMF, 
and  a  relationship  between  mass  and  main-sequence  lifetime  from  Adams 
(2008),  normalized  to  tms  =  10  Gyr  at  solar  mass.  To  be  consistent  with  the 
IMF  used  to  calibrate  the  eagle  simulations,  we  choose  the  minimum  and 
maximum  stellar  masses  to  be  Mmin  =  0.1  Mq  and  Mmax  =  100  Mq. 


Note  that,  since  the  time  at  which  observers  exist  is  irrelevant  to 
the  mass-weighted  measure,  the  ‘Star  formation  +  fixed  delay’  and 
‘Star  formation  +  main-sequence  lifetime’  models  give  identical 
results.  This  is  not  the  case  for  the  causal  patch  and  causal  diamond 
measures  -  a  later  observer  at  the  same  comoving  position  may  be 
outside  the  patch/diamond,  and  so  does  not  contribute  to  the  integral 
in  equation  (15). 

3.  Star  formation  +  metals:  the  raw  materials  for  life  are  the 
product  of  stellar-nucleosynthesis,  and  in  particular  metals  that  have 
been  ejected  from  stars  and  returned  to  the  interstellar  medium. 
Planets,  it  is  believed,  form  from  the  debris  discs  around  newly 
formed  stars,  and  stars  with  higher  metallicity  are  known  to  be 
more  likely  to  have  giant  planets  (Gonzalez  1997;  Fischer  &  Valenti 
2005).  However,  this  result  is  less  clear  for  smaller  rocky  planets 
(Buchhave  &  Latham  2015;  Wang  &  Fischer  2015).  There  must, 
of  course,  be  some  metallicity  dependence,  since  the  probability 
of  a  rocky  planet  forming  in  a  zero-metallicity  debris  disc  is  zero. 
We  make  the  simple  assumption  that  the  probability  of  a  rocky 
planet  forming  around  a  star  is  proportional  to  the  metallicity  of 
the  star-forming  gas,  Zsf,  so  that  the  observer  creation  rate  at  time 
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Figure  11.  The  ‘relative  observer  creation  rate’,  which  is  the  integrand  of  equation  (15)  or  (16).  The  value  of  A  is  shown  in  the  legend  in  the  top  left  panel. 
The  first  row  shows  the  mass-weighted  measure,  derived  directly  from  the  star  formation  rate  efficiencies  in  Fig.  5  and  metal  fraction  in  Fig.  6.  The  second 
row  shows  the  causal  patch  measure,  and  the  third  row  shows  the  causal  diamond  measure.  The  columns  show  the  different  observer  models.  The  first  column 
shows  the  star  formation  +  fixed  delay  model,  the  second  column  shows  the  star  formation  +  main-sequence  lifetime  model,  and  the  third  column  shows  the 
star  formation  -f  metals  model. 


t  is  proportional  to  the  number  of  planets  that  exist  around  main- 
sequence  stars, 

-^(Ooc  /  Zsf(0  Atar(^')/msa  -  (22) 

d^dV  Jo 

where  Zsf(^0  is  the  metallicity  of  star-forming  gas  at  time  t . 

5.3  Extrapolation 

The  integral  in  equations  (15)  and  (16)  is  over  all  of  cosmic  time, 
but  our  simulations  only  extend  to  a  finite  time.  They  capture  the 
initial  burst  of  star  formation  in  the  universe,  and  so  are  converging 
thanks  to  the  isolation  of  haloes  by  the  acceleration  of  the  expansion 
of  space.  There  will,  however,  be  a  trickle  of  star  formation  into  the 
future  in  our  galaxies,  which  our  simulations  do  not  capture.  Look¬ 
ing  at  the  decline  of  star  formation  in  the  A  >  Aq  x  10  simulations, 
we  extrapolate  our  simulations  using  an  exponential  decrease  in  star 
formation  (rate)  efficiency  (SFE)  with  time  [SFE  =  aexp(  —  bt)], 
for  constants  a  and  b  that  are  derived  from  the  final  few  Gyr 
of  the  simulation.  Eor  our  simulations,  Zsf(0  has  converged; 
extrapolating  by  fitting  an  exponential  makes  only  a  negligible 
difference. 


We  also  use  the  A  =  0  simulation  to  calculate  the  relevant  quan¬ 
tities  for  0  <  A  <  Aq.  Eor  A  <  Ao  X  0.1,  the  time  at  which  the 
universe  begins  to  accelerate  is  greater  than  the  limits  of  our  simu¬ 
lation,  at  which  time  the  star  formation  rate  efficiency  has  peaked 
and  is  declining.  By  using  the  A  =  0  simulation,  there  will  be  no 
difference  in  the  observer  model,  but  there  will  be  a  difference  in 
the  causal  patch  and  causal  diamond  because  these  depend  on  A. 
Note  that  the  both  of  these  measures  diverge  for  A  =  0,  so  we  only 
consider  universes  with  A  >  0. 

As  noted  in  Section  2.4,  in  the  far  future  of  our  simulations,  a 
variety  of  very  slow  (on  Gyr  time-scales)  processes  may  become 
relevant  but  are  not  captured  by  our  simulations.  In  particular,  in  the 
distant  future  of  the  A  =  0  cosmology,  the  growth  function  grows 
without  limit  and  it  is  likely  that  close  to  100  per  cent  of  the  mass 
in  the  universe  is  found  in  galaxies.  However,  in  an  old,  extremely 
diffuse  galactic  disc,  baryons  may  be  more  likely  to  accrete  directly 
on  to  dead  stellar  remnants  and  BHs  than  collapse  into  a  fresh 
star.  What  is  the  long-term  fate  of  the  interstellar  medium  in  our 
isolated  galaxies?  Further  modelling  may  be  able  to  derive  the 
expected  fraction  of  baryons  that  will  form  stars  into  the  distant 
future,  and  in  particular,  stars  that  are  likely  to  host  planets.  Here, 
in  the  absence  of  such  a  model,  we  will  simply  extrapolate  the 
simulations. 
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Figure  12.  The  relative  probabilities  per  unit  log  A  from  equation  (15)  or  (16);  each  line  integrates  a  panel  of  Fig.  1 1  over  cosmic  time.  The  left-hmd  panel 
shows  the  mass-weighted  measure,  the  middle  panel  shows  the  causal  patch  measure,  and  the  ng/zt-hand  panel  shows  the  causal  diamond  measure. 


Table  3.  Median  and  ‘In’  (68  per  cent)  probability  limits  of  the  cosmological  constant  for  the  three 
multiverse  measures  and  three  observer  models.  For  the  causal  patch  and  causal  diamond  measures,  the 
value  in  brackets  shows  the  median  value  of  the  cosmological  constant  using  the  observer  creation  rate 
(per  unit  mass)  from  the  A  =  0  simulation.  This  illustrates  the  effect  of  these  measures. 
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5.4  Predicting  the  cosmological  constant 

Fig.  11  shows  the  integrand  in  the  numerator  of  equation  (15)  or 
(16),  which  combines  the  observer  creation  rate,  the  measure,  and 
the  chosen  comoving  volume  (if  needed).  We  will  call  this  the 
‘relative  observer  creation  rate’.  The  value  of  A  is  shown  in  the 
legend  in  the  top  left  panel.  The  first  row  shows  the  mass-weighted 
measure,  derived  directly  from  the  star  formation  rate  efficiencies 
in  Fig.  5  and  metal  fraction  in  Fig.  6.  The  second  row  shows  the 
causal  patch  measure,  and  the  third  row  shows  the  causal  diamond 
measure. 

The  columns  show  the  different  observer  models.  The  first  col¬ 
umn  shows  the  star  formation  +  fixed  delay  model,  the  second 
column  shows  the  star  formation  +  main- sequence  lifetime  model, 
and  the  third  column  shows  the  star  formation  +  metals  model. 

The  first  row  (mass  weighted)  shows  most  directly  the  effects 
of  the  different  observer  models.  The  grey  Aq  x  0.01  model  is 
indistinguishable  from  the  Aq  x  0. 1  model  because  they  both  use 
the  results  of  the  A  =  0  simulation,  as  described  above.  In  the  second 
column,  we  can  see  the  effect  of  folding  in  the  main- sequence  stellar 
lifetime.  The  decline  in  the  observer  creation  rate  follows  the  decline 
in  the  main- sequence  fraction,  as  the  initial  burst  of  stars  formed 
in  the  first  10  Gyr  after  the  big  bang  grow  old.  Because  of  the 
abundance  of  small,  long-lived  stars,  observers  are  created  even  at 
very  late  times  in  the  universe  (Loeb,  Batista  &  Sloan  2016).  The 
addition  of  metal- weighting  further  diminishes  observer  creation  in 
large  A  universes,  as  they  have  fewer  stars  and  fewer  bound  metals 
to  make  planets  around  their  stars. 

The  second  and  third  rows  show  the  causal  patch  and  causal 
diamond  measures.  Note  that  the  Aq  x  0.01  and  Aq  x  0.1  curves  are 
distinguishable  because  of  the  difference  in  the  comoving  volume 
V{t\  A).  For  each  observer  model,  these  measures  show  similar 
trends.  The  smaller  comoving  volume  at  earlier  times  in  the  causal 
diamond  moves  the  peak  to  slightly  later  times,  but  otherwise  the 


two  measures  are  very  similar.  The  main  effect  of  these  measures 
is  to  decrease  the  relative  observer  creation  rate  exponentially  once 
the  expansion  of  the  universe  begins  to  accelerate.  This  somewhat 
cancels  out  the  effect  of  the  longer  main- sequence  lifetimes  in  the 
second  and  third  observer  models. 

Fig.  12  shows  the  relative  probabilities  p{A\MB)  from  equation 
(15)  or  (16);  each  line  integrates  a  panel  of  Fig.  1 1  over  cosmic  time. 
We  plot  the  probability  per  unit  log  A,  and  normalize  by  setting  the 
maximum  value  to  one,  rather  than  integrating  over  the  limited 
range  of  A.  Note  that  in  the  left  plot  (mass- weighted  measure),  the 
‘SF  +  delay’  and  ‘SF  +  lifetime’  curves  are  indistinguishable  -  the 
integral  over  cosmic  time  is  not  affected  by  the  5  Gyr  delay,  and 
cancels  out  the  effect  of /ms-  For  the  range  of  A  we  consider,  the 
median  and  ‘la’  (68  per  cent)  values  are  shown  in  Table  3. 

As  we  have  noted  previously,  the  decline  in  star  formation  in  our 
universe  after  t  =  3.5  Gyr  is  not  due  to  the  effect  of  the  cosmologi¬ 
cal  constant  (Salcido  et  al.  2018).  Universes  without  a  cosmological 
constant  show  a  similar  decline.  The  initial  burst  of  star  formation 
in  the  universe,  then,  is  not  dramatically  affected  by  moderate  in¬ 
creases  in  A.  Only  for  A  >  Aq  x  30  do  we  see  a  significant  effect  on 
the  total  number  of  stars  in  the  universe.  Thus,  in  the  mass-weighted 
measure,  the  probability  distribution  for  A  is  reasonably  fiat  to  large 
values  ('^Aq  x  30).  The  median  value  in  this  case  is  60  times  larger 
than  the  observed  value.  Adding  metal-weighting  to  the  observer 
model  increases  the  suppressing  effects  of  A,  but  the  median  value 
is  still  ~45  times  larger  than  the  observed  value.  While  these  dis¬ 
tributions  are  broad,  most  of  the  probability  is  at  large  values  of  A. 
Table  4  shows  the  probability  that  the  cosmological  constant  ob¬ 
served  by  a  typical  observer  is  less  than  or  equal  to  the  value  in  our 
Universe  (Aq)  for  the  three  multiverse  measures  and  three  observer 
models.  For  the  mass- weighted  measure,  this  probability  is  small 
(2  per  cent). 

Note  that  the  results  above  are  for  the  parameter  range 
Aq  X  0.01  <  A  <  Aq  X  300.  The  results  for  small  values  of 
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Table  4.  The  probability  that  the  cosmological  constant  observed  by  a 
typical  observer  is  less  than  or  equal  to  the  value  in  our  Universe  (Aq)  for 
the  three  multiverse  measures  and  three  observer  models.  For  the  causal 
patch  and  causal  diamond  measures,  these  probabilities  are  greater  than 
50  per  cent,  but  the  value  for  the  mass-weighted  measure  is  small. 


Prob  A  <  Aq 

Mass  weighted 

Causal  patch 

Causal  diamond 

SF  +  delay 

1 .9  per  cent 

86  per  cent 

73  per  cent 

SF  +  lifetime 

1 .9  per  cent 

90  per  cent 

86  per  cent 

SF  +  metals 

2.5  percent 

93  per  cent 

90  per  cent 

A  have  converged,  but  increasing  the  upper  limit  increases  the  me¬ 
dian  value.  If  we  extrapolate  the  probability  distribution  to  larger 
values  of  A,  we  find  that  the  median  value  for  A  is  ~200  Aq  for  the 
mass-weighted  measure. 

For  the  causal  patch  and  causal  diamond  measures,  the  fact  that 
the  comoving  volume  in  the  measure  decreases  with  time  sup¬ 
presses  large  values  of  the  cosmological  constant,  independently 
of  the  effect  on  the  observer  creation  rate  d^Ttobs/dMV.  This  leads 
to  the  >  50  per  cent  probabilities  for  small  values  of  the  cosmo¬ 
logical  constant  (Table  4).  To  illustrate  the  effect  of  the  measure, 
we  calculate  the  median  value  of  the  cosmological  constant  using 
the  observer  creation  rate  (per  unit  mass)  from  the  A  =  0  simu¬ 
lation.  The  results  are  shown  in  brackets  in  Table  3.  These  values 
are  consistently  larger  than  the  actual  median  value,  but  only  by  a 
small  factor.  Thus,  the  causal  patch  and  causal  diamond  measures 
are  playing  the  dominant  role  in  setting  the  expected  value  of  the 
cosmological  constant.  The  predicted  value  of  A  is  set  by  the  time 
at  which  the  star  formation  efficiency  peaks  in  universes  with  small 
values  of  A,  which  is  set  by  other  cosmological  and  physical  pa¬ 
rameters.  The  decline  of  star  formation  efficiency  with  A  plays  a 
secondary  role. 

Our  results  are  broadly  consistent  with  the  analytic  model  of 
Bousso  &  Leichenauer  (2010),  who  find  that  for  fixed  values  of  the 
primordial  inhomogeneity  Q  and  spatial  curvature,  and  for  A  >  0, 
the  causal  patch  and  causal  diamond  measures  predict  a  value  of 
0.1  <  A/Aq  <  10,  depending  on  the  model  for  observers.  As 
noted  there,  the  suppression  of  structure  formation  by  accelerating 
expansion  is  only  important  for  a  cosmological  constant  of  the  of 
order  Aq  x  100.  Thus,  the  agreement  between  our  calculations 
is  due  to  the  ‘geometric’  effects  of  the  causal  patch  and  causal 
diamond  measures;  the  astrophysics  of  galaxy  formation  does  not 
prefer  values  of  the  cosmological  constant  less  than  Aq  x  100. 

6  CONCLUSIONS 

Models  of  the  very  early  universe,  including  inflationary  models,  are 
argued  to  produce  varying  universe  domains  with  different  values 
of  fundamental  constants  and  cosmic  parameters.  In  such  models, 
predicting  observations  necessarily  involves  understanding  where 
observers  are  created  in  the  multiverse.  In  particular,  this  anthropic 
approach  has  been  used  to  predict  the  value  of  the  cosmological 
constant. 

Using  the  cosmological  hydrodynamical  simulation  code  from 
the  EAGLE  collaboration,  we  have  investigated  the  effect  of  the  cos¬ 
mological  constant  on  the  formation  of  galaxies  and  stars.  This 
SPH  code  follows  the  gravitational  collapse  of  matter  in  an  expand¬ 
ing  universe,  incorporating  sub-grid  recipes  for  radiative  cooling 
for  1 1  elements,  star  formation,  stellar  mass-loss,  energy  feedback 
from  star  formation,  gas  accretion  on  to  and  mergers  of  supermas- 
sive  BHs,  and  AGN  feedback.  We  simulate  universes  with  values 
of  the  cosmological  constant  ranging  from  A  =  0  to  Aq  x  300, 


where  Aq  is  the  value  of  the  cosmological  constant  in  our  Uni¬ 
verse.  For  larger  values  of  the  cosmological  constant,  the  time  at 
which  the  expansion  of  the  universe  begins  to  accelerate  declines  as 
fA  oc  (A/Ao)”'/^. 

Our  Universe  shows  a  peak  in  the  global  star  formation  rate  at 
t  =  3.5  Gyr,  coming  after  the  peak  in  the  halo  matter  accretion 
rate  at  t  =  I  Gyr.  By  the  time  the  expansion  of  our  Universe  be¬ 
gins  to  accelerate  (at  t  =  1.6  Gyr),  the  global  halo  mass  accretion 
rate  has  dropped  to  about  10  per  cent  of  its  earlier  maximum,  and 
most  of  the  mass  that  will  ever  accrete  into  haloes  has  already 
accreted.  As  a  result,  increases  in  A  of  even  an  order  of  magni¬ 
tude  have  a  small  effect  on  the  star  formation  efficiency  of  the 
universe. 

One  interesting  effect  that  affects  the  raw  materials  of  life  is 
stellar  and  AGN  feedback.  In  our  Universe,  these  processes  slow 
star  formation  by  sending  baryons  back  into  the  outer  parts  of  the 
halo  and  the  local  intergalactic  medium.  This  material  is  largely 
recycled  into  the  galaxy  after  ~1  Gyr,  and  forms  a  later  generation 
of  stars.  But  in  universes  with  A  >  Aq  x  10,  much  of  this  material 
is  lost  to  the  intergalactic  medium,  carried  away  by  the  accelerating 
expansion  of  the  universe  rather  than  reaccreting.  The  net  baryon 
accretion  rate  becomes  negative  as  more  material  is  lost  to  galactic 
winds  than  is  accreted/reaccreted. 

In  universes  with  larger  values  of  A,  galaxies  quickly  become 
isolated  from  their  cosmic  surroundings.  The  familiar  ecosystem  of 
galaxies  in  our  Universe,  which  balance  accretion,  major  and  minor 
merging,  galactic  cannibalism,  star  formation,  galactic  winds,  and 
reaccretion,  is  reduced  to  a  closed  box,  as  galaxies  become  island 
universes,  surrounded  by  vacuum  and  isolated  from  the  rest  of  the 
matter  in  the  universes.  They  burn  through  their  finite  matter  supply, 
forming  stars  at  a  decreasing  rate. 

We  use  our  simulations  to  predict  the  observed  value  of  the  cos¬ 
mological  constant,  given  a  measure  of  the  multiverse.  We  consid¬ 
ered  three  simple  but  plausible  models  for  where  we  would  expect 
observers  to  be  created  in  our  simulations,  and  three  measures  of 
the  multiverse. 

In  the  mass-weighted  measure,  with  a  uniform  probability  that 
a  given  mass  element  in  the  universe  will  inhabit  a  region  with 
a  given  value  of  the  cosmological  constant,  the  predicted  size  of 
A  is  determined  by  the  decline  in  the  star  formation  efficiency  of 
the  universe.  For  the  reasons  described  above,  this  is  relatively  flat 
as  a  function  of  A,  and  so  the  predicted  (median)  value  is  50- 
60  times  larger  than  the  observed  value.  The  probability  of  observ¬ 
ing  a  value  as  small  as  our  cosmological  constant  Aq  is  ~2  per  cent. 
In  this  case,  an  anthropic  argument  for  value  of  A,  while  doing 
much  better  than  the  famous  120  orders-of-magnitude  discrep¬ 
ancy  from  quantum  field  theory,  is  not  a  particularly  successful 
prediction. 

For  the  causal  patch  and  causal  diamond  measures,  which  con¬ 
sider  a  subset  of  the  universe  that  depends  on  A,  the  predicted  value 
is  within  a  factor  of  a  few  of  the  observed  value.  But  this  has  very 
little  to  do  with  the  decline  in  the  star  formation  efficiency  (and  so, 
presumably,  observer  creation  rate)  with  A.  It  is  a  result  of  the  rapid 
decrease  in  the  size  of  the  causal  patch/diamond  with  increasing 
cosmological  constant. 

We  stress  again:  this  is  no  reason  to  prefer  the  causal  patch  and 
causal  diamond  measures.  This  is  not  an  observational  test  of  these 
measures.  A  specific  multiverse  model  must  justify  its  measure  on 
its  own  terms,  since  the  freedom  to  choose  a  measure  is  simply  the 
freedom  to  choose  predictions  ad  hoc. 

We  conclude  that  the  impact  of  the  cosmological  constant  on  the 
formation  of  structure  in  the  universe  does  not  straightforwardly 
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explain  the  small  observed  value  of  A.  The  prediction  depends 
crucially  on  the  measure.  If  the  observer  creation  rate  had  been 
sufficiently  sharply  peaked  at  values  near  Aq,  the  measure  would 
not  much  matter.  But  in  fact,  in  the  absence  of  a  multiverse  model 
that  can  convincingly  justify  a  measure,  it  is  not  clear  whether  the 
anthropic  prediction  A  is  successful.  Future  work  will  consider 
varying  more  cosmological  and  fundamental  parameters  to  shed 
more  light  on  which  kind  of  universe  is  to  be  expected  from  a 
multiverse. 
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